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Abstract 

Considering interacting (antiferromagnetically correlated) electrons, we de- 
rive a generalized Harper's equation for the square lattice of infinite size. We 
obtain an analytic expression for the density of states from the newly derived 
Harper's equation. We present a predicted phase diagram of staggered magne- 
tization in the plane of temperature vs doping rate and discover a possibility 
of reentrant behavior of the staggered magnetization even in the presence of 
applied magnetic field. It is shown that below a critical electron correlation 
strength (Coulomb repulsion) the staggered magnetization in the presence of 
magnetic field vanishes at an even denominator q value but not at odd q of a 
given magnetic flux quantum per plaquette, p/q. 

PACS numbers: 75.30.Kz, 75.50.Ee, 71.10.Fd 
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I. INTRODUCTION 



Since the discovery of high temperature superconductors and the related insulating ma- 
terials, there has been a lot of interest in two-dimensional spin-| magnetic or antiferromag- 
netically correlated electron systems.0! However not much attention has been paid to the 
magnetic properties of these systems coupled to an external magnetic fieldJHl The original 
Harper's equation@~! is concerned with the energy dispersion involving the systems of non- 
interacting electrons under the magnetic field. Hence it is of great interest to study how the 
systems of interacting (antiferromagnetically correlated) electrons behave under an external 
magnetic field. We derive a generalized Harper's equation which describes the dispersion 
of antiferromagnetically correlated electrons under the applied magnetic field. Earlier we 
paid attention to the dispersion of the antiferromagnetically correlated electrons only at half 
filling (and thus with no doping) and at zero temperature, by considering the square lattice 
of finite size.i On the other hand, in the present study we derive a generalized Harper's 
equation for the square lattice of the infinite size, and examine the dispersion relation as a 
function of temperature and doping rate. An analytic expression for the density of states 
is obtained from the generalized Harper's equation for the system of antiferromagnetically 
correlated electrons. We present a phase diagram of staggered magnetization in the plane 
of temperature vs doping rate and find the hitherto-unnoticed reentrant behavior of the 
staggered magnetization even in the presence of external magnetic field. Finally it is shown 
that below a critical electron correlation strength the staggered magnetization disappears 
at an even denominator value of q but not at odd q of a given magnetic flux quantum per 
plaquette, p/q. 

II. GENERALIZED HARPER'S EQUATION AND DENSITY OF STATES 

We write the Hubbard model Hamiltonian describing the two-dimensional system of 
antiferromagnetically correlated electrons under an external magnetic field,i 
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where t is the hopping integral; A, the electromagnetic vector potential; <po = — , the 
elementary flux quantum; Z7, the on-site Coulomb repulsion energy, and /i, the chemical 
potential, (ij) stands for summation over nearest neighbor sites i and j. c\ a (c ia ) is the 
creation (annihilation) operator of an electron of spin a at site i, and (n^), the number 
operator of an up-spin (down-spin) electron at site i. 

The staggered magnetization (antiferromagnetic order) at site i is written as rrij = 
giQ-ri a(cl a c ia ), where Q = (n,ir) and = (i x ,i y ) with i x and i y being integers with 
the lattice spacing of unity. Introducing a uniform staggered magnetization m and a uni- 
form doping rate 5, i.e., 

™ = ^E etQ 'M40 , (2a) 

^£k> , (26) 

with the number of lattice sites N, and using the Landau gauge A = B(0, x, 0), we obtain 
the mean field (Hartree-Fock) Hamiltonian in the momentum space, 

H = - t H [ 2 cos k x cl a c ka + e- ik »cl_ gja c ka + ^ c£ +g)(J c k(7 

kcr 

-^E aC k+Q, CT C k(7 + [y(l -5) -fi]J2 C Lcka , (3) 
ktr ktr 

where g = \ 2ir4^, o) = (27r£, with |, the number of flux quanta per plaquette. The first 
bracketed term in (0) represents hopping processes; the first term in the bracket represents 
the nearest neighbor hopping in the x-direction and the last two terms in the bracket, the 
nearest neighbor hopping in the ^-direction. Because of the choice of the Landau gauge 
A = B(0,x,0), the electron acquires no phase when it hops in the x-direction, while it 
acquires a phase when it hops in the ^-direction, and the electromagnetic vector potential 
A shifts the wave vector of electron in the fc^-direction by g = |g| = 27T-. The second term 
results from the nature of the antiferromagnetic spin order of correlated electrons, which 
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causes the wave vector to shift by Q. The last term represents energy shift by ~(1 — 5) as 
a result of hole doping. 

The Hamiltonian (EJ) can be written as 



H = H + H 1 



r U 



kcr 



where 
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with M n = 2cos(/c x + ng). The summation J2'is. * s over the reduced Brillouin zone, 
{(k x ,k y )\ — - < k x < -, — ~ < k y < -}. The Brillouin zone is reduced by 1/q be- 
cause there exist q plaquettes per magnetic unit cell, and is further reduced by 1/2 as a 
result of staggered magnetization (or antiferromagnetic order). The diagonal matrix Tj. 
is associated with electron hopping and contains information on the phase modulation of 
hopping electrons under the influence of the external field. The off-diagonal matrix V CT 
represents the antiferromagnetic electron correlation. 

From the eigenvalue equation of the Hamiltonian matrix H^o- with the identity matrix I, 

det(H fcCT - E k T) = , (5) 

we obtain the quasiparticle energy dispersion Ef. of antiferromagnetically correlated electrons 
in the presence of magnetic field. In the limiting case of noninteracting electrons (U = 0) 
the 'generalized' Harper's equation (||) above is reduced to the original Harper's equation 
derived by Hasegawa et al.i, that is, 

det(T fc - e fc I) =0 , (6) 

where Ek is the energy dispersion of noninteracting electrons in the presence of magnetic 
field. Following Hasegawa et al.i, Eq. @ can be rewritten in a further simplified form, 

7(e) = cos(qk x ) + cos(qk y ) , (7) 

where 7(e) is given in Table | for various values of p/q including Hasegawa et al.'s results.^ 
Now for the case of antiferromagnetically correlated electrons we obtain from the diagonal- 
ization of the Hamiltonian matrix in Eq. (4c) above, 
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with the band gap, 2A = mU. The band gap 2A is seen to depend on the magnitude of both 
the staggered magnetization m and the electron correlation strength (Coulomb repulsion) 
U. 

The energy dispersion relation Eq. (||) above leads to the density of states, 
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where |e| = y E 2 — ( !2 f^) and K is the complete elliptic integral of the first kind,0 

K(a) 



y 1 — a 2 sin 2 

Eq. @ above represents the density of states for the systems of correlated electrons in the 
presence of magnetic field. In the limiting case of vanishing electron correlation, i.e., {7 = 0, 
it becomes exact. In Fig. [I] we display the density of states with {7 = predicted from the 
analytic expression of the density of states in Eq. @. Encouragingly the analytic result 
is in excellent agreement with the numerical results obtained by Hasegawa et al.i In the 
presence of magnetic field with the flux quanta per plaquette p/q, a single band splits into 
q subbands. This is analogous to the energy level splitting into Landau levels for electrons 
embedded in a continuum state under a magnetic field. 

III. PHASE DIAGRAM AND REENTRANT BEHAVIOR OF STAGGERED 
MAGNETIZATION IN APPLIED MAGNETIC FIELDS 



The staggered magnetization m and the chemical potential fi vary with both the tem- 
perature T and the doping rate 5; they are obtained from the following self-consistent mean 
field equations, which are derived from the use of Eqs. (2), 
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where £ (±) = ±£ + f (1 - 5) - //. 

We examine the dependence of staggered magnetization on the external magnetic field 
in the T-8 plane,0 by numerically solving Eqs. (10) for m and fi at each temperature T 
and doping rate 5. The phase diagram of the staggered magnetization in the T-8 plane is 
displayed in Fig. |2| for several values of p/q. Interestingly we find that the reentrant behav- 
ior of the staggered magnetization appears even in the presence of the external magnetic 
field. However, the reentrant behavior of staggered magnetization in the absence of mag- 
netic field has earlier been discovered by other investigators.0~@ Halvorsen et al.0 found the 
reentrant behavior using the Hubbard model within the self-consistent second-order weak 
[/-perturbation treatment. Their studies are limited to a narrow range with small U com- 
pared to our present approach which can deal with the entire range of U. Recently Inaba 
et alH found a similar reentrant behavior using the t-J Hamiltonian in the slave-boson 
representation. However their studies refer to the case of a large U limit due to the use of 
the t-J Hamiltonian. Unlike our present study their studies above refer to the reentrant 
behavior in the absence of the external magnetic field. 

We now investigate the reentrant behavior of the staggered magnetization in detail. In 
Fig. [3] we display the staggered magnetization m at p/q = 1/2 as a function of temperature 
T at various doping rates 5. At half filling (5 = 0) the staggered magnetization reaches a 
maximum value at zero temperature. On the other hand, away from half filling (5 ^ 0) the 
predicted staggered magnetization shows a maximum at a finite temperature. Above this 
temperature the reentrant behavior of a paramagnetic state is predicted. The cause of the 
reentrant behavior can be explained from the nesting property of the energy surface at saddle 
points.lli0 In Fig. £| we display the variation of energy dispersion of the highest occupied 
subband with the external magnetic field (or p/q). The saddle points in the Brillouin zone 
are denoted by black dots at the bottom of the graph. The adjacent saddle points are 
separated by the well-defined nesting vector of Q/q, by which the staggered magnetization 
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is defined. In Fig. [5] we display the variation of Fermi surface (thick solid lines) with the 
magnetic field at finite doping rates. As the temperature increases, the Fermi surface tends 
to smear out, which causes increment in the number of nesting channels, and consequently 
the staggered magnetization arises. As the temperature still increases, further smearing of 
the Fermi surface opens other channels than the nesting channels. As a result the staggered 
magnetization will eventually disappear to allow transition to a paramagnetic state. This 
feature is well depicted in Figs. || and |[ 

IV. STAGGERED MAGNETIZATION AT ZERO TEMPERATURE AND AT 
HALF FILLING IN A MAGNETIC FIELD 

Now we investigate the staggered magnetization at half filling 5 = and at T = OK in a 
magnetic field. The chemical potential is given by \i = ^, and Eq. (10b) is trivially identified 
since the right hand side of the equation becomes zero, thus satisfying the condition of half 
filling, that is, 5 = as it should be. Eq. (10a) can be expressed as 

* 1 (ii) 

Here go(e) is the density of states of noninteracting electrons in the presence of magnetic 
field. It is easily obtained from Eq. @. Eq. (|TT|) is in a similar form to the gap equation 
that appears in the spin density wave theory of cuprate materials. We will use it for the 
determination of the staggered magnetization m for the system of interacting electrons (with 
correlation strength U) in the presence of magnetic field. In Fig. ^ the oscillatory staggered 
magnetization is displayed as a function of magnetic field (specifically p/q), for several chosen 
values of correlation strength U. The solid lines are the results of self-consistent calculations 
for a 20 x 20 finite square lattice.i Various other symbols represent the results from the 
newly derived analytic relation (|IT]) above. Encouragingly they are in good agreement with 
the self-consistent calculations for the finite size square lattice. 

In the following we explain the oscillatory behavior of staggered magnetization observed 



in our earlier workrl At even denominator values of q in p/q the staggered magnetization 
is predicted to disappear (e.g., see the case of p/q — 1/2). This feature is well depicted in 
Fig. |6|. We now define the critical electron correlation strength (Coulomb repulsion) U p / q 
as a value below which the staggered magnetization vanishes, i.e., m = 0. The predicted 
staggered magnetization from Eq. flTT| ) vanishes at the even q values below a critical value 
Up/g, i.e., U < Up/ q . On the other hand in the absence of magnetic field the staggered 
magnetization tends to appear even at small values of £7 ill The critical correlation strength 
Up/q is obtained by substituting m = in Eq. ([□]), and is shown for various values of p/q in 
Table [II]. F° r °dd q in p/q, the integral in Eq. (|ll|) is logarithmically divergent, and thus U p / q 
does not exist. Although not numerically precise for the case of finite size calculations (solid 
lines) J§ with the even denominator values of q in p/q a propensity of vanishing staggered 
magnetization is seen below U p / q as shown in Fig. || Thus the oscillatory behavior of the 
staggered magnetization is found to occur owing to its disappearance distinctively at the 
even denominator values of q in p/q. 



V. CONCLUSION 

We derived a generalized Harper's equation for the energy dispersion relation of inter- 
acting (antiferromagnetically correlated) electrons in an external magnetic field. Unlike 
the original Harper's equation which deals only with noninteracting electron systems, the 
generalized Harper's equation derived in Eq. (|5]) (with Eq. (4c)) has a definite merit of 
studying the physical properties of correlated electron systems in the presence of the ex- 
ternal magnetic field. From this Harper's equation we derived an analytic formula for the 
density of states of the antiferromagnetically correlated electrons in the magnetic field. For 
the limiting case of noninteracting electrons, the analytic equation for the density of states 
is found to be in good agreement with the numerical work of Hasegawa et al.i Further we 
presented the phase diagram of staggered magnetization in the plane of temperature T vs 
doping rate 5 as a function of magnetic field (specifically, flux quanta per plaquette p/q) 
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and correlation strength U. From this study we demonstrated a possibility of reentrant 
behavior of staggered magnetization even in the presence of the applied magnetic field. A 
more accurate account of electron correlations beyond the mean field approximations may 
not alter the qualitative nature of the reentrant behavior that we discovered in this study. 
Although not reported here, we find from the exact diagonalization study of Hubbard model 
that the accurate account of correlations does not affect the qualitative finding here. In 
the present study we neglected the Zeeman coupling. At such large Coulomb repulsion en- 
ergies as U = 8t which we used in our calculations, we find that the Zeeman effect does 
not substantially alter the observed staggered magnetization. For the case of half filling 
at zero temperature, we obtained a gap equation (|TT| ) for the determination of staggered 
magnetization for antiferromagnetically correlated electron systems at a given correlation 
strength U in the presence of magnetic field. From this derivation we were able to deter- 
mine a critical correlation strength U p / q below which staggered magnetization disappears at 
even denominator values of q (but not at odd denominator values of q) of the magnetic flux 
quanta per plaquette, p/q. 
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TABLE CAPTIONS 



TABLE I. 7(e) for various values of magnetic flux quanta per plaquette p/q. The energy dis- 
persion of noninteracting electrons at a given p/q is determined from the Harper's 
equation, 7(e) = cos(qk x + qk y ). 

TABLE II. Critical correlation strength U p / q as a function of magnetic flux quanta per plaquette 
p/q. 
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FIGURE CAPTIONS 



FIG. 1. Density of states of noninteracting electrons in the presence of applied magnetic field 
for various flux quanta per plaquette p/q based on the analytic expression @. 

FIG. 2. Phase diagram of staggered magnetization for several values of magnetic flux quanta 
per plaquette p/q. Each line represents a boundary between the antiferromagnetic 
(AF) phase and the paramagnetic (PM) phase. 

FIG. 3. Variation of staggered magnetization as a function of temperature for various doping 
rates at p/q = 1/2 and U = 8t. 

FIG. 4. Variation of quasiparticle energy dispersion surface of the highest occupied subband 
with p/q. Black dots denote the saddle points of the surface. Rectangles at the 
bottom of the graph represent the reduced Brillouin zones, and arrows, the nesting 
vectors between two adjacent saddle points. 

FIG. 5. Variation of Fermi surfaces with p/q at zero temperature and at hole doping rates, (a) 
5 = 0.33 (b) 5 = 0.30 (c) 5 = 0.10 (d) 5 = 0.17. Black dots denote the saddle points; 
rectangles, the reduced Brillouin zones, and arrows, the nesting vectors between two 
adjacent saddle points. 

FIG. 6. Staggered magnetization (antiferromagnetic order) at T = OK as a function of p/q, 
the magnetic flux quanta per plaquette. The solid lines denote the results from self- 
consistent calculations on a 20 x 20 finite square lattice, and various other symbols 
represent the results calculated from the analytic expression ( |TT|) corresponding to the 
square lattice of infinite size. 
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TABLES 



27(f) 

e 2 -4 
—e 3 + 6e 
e 4 - 8e 2 + 4 
-e 5 + 10e 3 + e (-15 + 10 cos(^)) 
-e 5 + 10e 3 + e (-20 - 10 cos(f )) 
e 6 - 12e 4 + 245 2 - 4 
-e 1 + 14e 5 + (-49 + 14 cos(^)) e 3 + (42 - 28 cos(^) - 28 cos(^)) e 
-e 1 + 14e 5 + (-49 - 14 cos(^)) e 3 + (28 - 28 cos(^) + 56 cos(f )) e 
-e 7 + 14e 5 + (-56 - 14 cos(f ) + 14 cos(f )) £ 3 + (56 + 56 cos(^) - 28 cos(f )) e 
e 8 - 16e 6 + (72 - %y/2) e 4 + (-96 + 32^) e 2 + 4 
£ 8 - 16e 6 + (72 + 8v^) e 4 + (-96 - 32^) e 2 + 4 
-£ 9 + 18^ 7 + (-99 + 4v/3 cos(^) + 6 cos(^) + 4^008^1)) e 5 
+ (l86 - 24 v / 3cos(f|) -36 cos(f ) +4VScos(^) - 28 v / 3cos(g) + 24 cos(f )) £ 3 
+ (-126 + 36 v / 3cos(^) + 36cos(^) + 36^008(^1) - 54cos(f )) e 
-e 9 + 18e 7 + (-99 + 4-s/3cos(f|) - 4^008(^1) + 6 cos(f )) £ 5 
+ (l86-4 v / 3cos(^) -24 cos(f ) - 28>/3cos(f|) + 24^008(^1) -60 cos(f )) e 3 
+ (-126 + 54 cos(f ) + 36^ cos(ff) - 36^08(^1) + 90 cos(f )) e 
-e 9 + 18e 7 + (-99 - 4 v / 3cos(^) - 6 cos(f ) - 4 v / 3cos(f§) - 6 cos(^)) e 5 
+ (186 + 28 v / 3cos(^) + 60 cos(f ) + 24^008(11) + 4 v / 3cos(g) + 36 cos(f )) e 3 
+ (-126 - 36^ cos(^) - 90 cos(^) - 36^ cos(ff) - 36 cos( 
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TABLE I. 
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TABLE II. 
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